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In a previous paper we examined a geometric measure of entanglement based on the minimum 
distance between the entangled target state of interest and the space of unnormalized product 
states. Here we present a detailed study of this entanglement measure for n— qubit target states 
\ that are invariant under the permutation of any two qubits. We analytically obtain the permutation 

invariant unnormalized direct product states that extremize the distance function. We then solve 
{^JQ' for the Hessian to show that, up to the action of trivial symmetries, the solutions correspond to 

'~{ I local minima of the distance function. In addition, we show that the conditions that determine the 

, extremal solutions for general target states can be obtained directly by parametrizing the product 

in 

Ph. 



0^ 



states via their Schmidt decomposition. 
PACS numbers: 03.65.Ud, 03.67.Mn 



I. INTRODUCTION 



Given the vital role that quantum entanglement is thought to play in quantum computing Q 
and its role in quantum phase transitions[5], it is important to have at one's disposal a quantitative 
measure of entanglement. General conditions that must be satisfied by useful entanglement measures 
ps) _ were presented in In that paper a class of good measures was constructed based on the notion of 
^ ■ minimum distance between the state in question (i.e. the target state) and the nearest disentangled 
state. There has been much work on so-called geometric measures of entanglement based on this 
distance to the nearest product state [T-II . Geometrical entanglement was applied to multipartite 



^ (3 qubit) systems in [12i] and related to the notion of entanglement witness in (isl . There 
^ has more recently been a resurgence of interest in various aspects of geometrical entanglement 



including formalism, generalizations and calculations for specific systems (see |15l-|24|. for a partial 
list). Explicit calculations of such measures for multipartite systems are generally very complicated 
and difficult to compute, except for target states with a high degree of symmetry. 



^ ' In a previous paper |25|, we studied a variation of the geometric measure of entanglement based 
5-1 on the minimum distance between an unnormalised product state and a target entangled state |28|. 

- The equations determining the (unnormalized) geometrical entanglement of a multipartite state were 
presented for a general target state. As shown in detail in 125*1 and reviewed at the end of the next 
section, there is a simple geometrical relationship between the geometrical entanglement calculated 
using unnormalized product states and those that are normalized. Although the equations are non- 
linear and hence difficult to solve in general, for target states with a large degree of symmetry it is 
possible to find analytic solutions that extremize the distance function. 

For a given solution to the variational equations to yield a meaningful measure of entanglement, 
it must be a local minimum of the distance function[29], and not a local maximum, or an infiection 
point in any direction within the parameter space. This information is provided by the Hessian for 
the system, which is the matrix of second variational derivatives evaluated at the extremum. In 
particular, all nonzero eigenvalues of the Hessian must be positive. Moreover, any zero eigenvalues 
should correspond to trivial symmetries of the system. 
In this paper, we address this issue in the context of a maximally symmetric, permutation invariant 
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target state consisting of q qubits In addition, we are able able to show quite generally that the 
geometric measure of entanglement has a direct interpretation in terms of a Schmidt decomposition of 
the multi-partite product states that enter the distance function. This result extends the connection 
with the Schmidt decomposition that was established in [25] for bi-partite systems. 

The paper is organized as follows: in the next Section we establish notation and review the 
geometric entanglement measure and resultant extrema conditions introduced in [25i]. Section [XTTl 
presents a general discussion of the Hessian and evaluates it for a maximally symmetric target state of 
q qubits. Section HVl presents the eigenvalues and eigenvectors of the Hessian (the explicit derivation 
is delegated to Appendix A) while Section |V] discusses the physical/geometrical interpretation of the 
eigenvectors in terms of their action on the space of product states. In Section |VT] we use the results 
of the previous sections to derive the analytic solution for the minimum distance and verify explicitly 
that it is a local minimum. In Appendix B, we illustrate that while our symmetric ansatz for the 
product states does solve the extremization conditions for a target state with less symmetry, it does 
not in that case correspond to a local minimum. Section IVIII describes the our distance measure in 
terms of a parametrization using the Schmidt decomposition of the product states. We end with 
Conclusions and prospects for future work. 



II. NOTATION AND SUMMARY OF PREVIOUS RESULTS 

We start by defining our notation. We consider a system of q qubits. The dimension of the 
corresponding Hilbert space is n = 2''. We decompose the system into a set of q subsystems, each of 
dimension 2. The subsystems are labelled A, B,C, ■ ■ ■ An arbitrary set of basis states of system A is 
labelled the basis states of B are \j), the basis states of C are \k), etc, and {i,j, A;, ■ ■ ■} G {0, 1}. 
Using this notation we write: 

\A)=aS), \B)=b,\j), \C) = Ck\k), ... (1) 

where the summation convention is implied. The coefficients ai,bj, .. are in general complex but we 
will henceforth consider them to be real for simplicity. Much of the following can be generalized to 
complex coefficients in a straightforward manner. 
The wave-function of an arbitrary normalised entangled pure state is written: 

\i^) = X^Jk■■■\^)(S)\J)^\k)■■■ , (^|V') = 1- (2) 

We consider a general product state of the form: 

10) = \A) ® ® |C) o . . . = ai\i) ® bj\j) ® Ck\k) (g) . . . = (f)ijk-\i) ® \j) ® ■ • ■ (3) 

Note that the state |0) is not assumed to be normalised: 

{<f)\<j)) = NANBNc... , NA = {A\A)=a*a,, Nb = {B\B) = b*b, , ■■■ (4) 

The distance between the states lip) and |0) is (25| : 

D'^ = {ip- 01?/' - 0) = \aibjCk ... - Xijk...\^ ■ (5) 

This distance function depends on the 2q parameters {aj, 6^, Cfc ■ ■ ■}. Taking the first derivatives with 
respect to the coefficients and setting the result to zero produces 2 equations which can be written: 

Az^2 ^ 2 {a,b,Ck ... - x^,k..) bjCk . . . = 2 (a, 6^ ... - 2 (x.,fc...&,c,. . . .) = , 2 G {0, 1} . (6) 
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In the same way, we can take derivatives with respect to the variables bj, Ck, ■ ■ ■ and set the resulting 
expressions to zero. We obtain a set of 2q equations which depend on the 2q variables {oj, bj,Ck- ■ ■}■ 
Substituting these solutions back into the distance function gives the condition that must be satisfied 
by an extremal solution 25j: 



D 



NaNbN, 



c ■ 



(7) 



We can write this condition in terms of a critical angle 6c which is defined as the angle between 
and \(h) at the extrema: 



Di 



COS 9r 



cos^ 9r 



VNaNbNc 



c 



critical 



The relationship of the unnormalized measure to the measure using a normalized product state is 
revealed by noting as in 25| that the minimization of the distance used to obtain the latter can 
be expressed in terms of the same variational principle as the former, but with the normalization 
condition on the product states imposed using a la gran ge multiplier. The corresponding result for 
the distance to a closest normalized product state is|25|: 



W - ^|07V - ^) = Dl + (1 - ^/Wi? 



2Dl 



(9) 



The final relationship above follows from a simple geometrical identity. The expectation therefore 
is that the two measures would in most cases be physically equivalent. One advantage of the 
unnormalized measure is that the variational equations without the normalization contraint are in 
principle somewhat simpler (although still non-linear and difficult to solve). 



III. THE HESSIAN 



The Hessian is the 2q x 2q matrix of second derivatives of the distance function with respect to 
the parameters {oj, 6^, ■ ■ ■}, evaluated at the given extremum. Denote by the complete set 
of numbers {a', 6^,...} parametrizing the unnormalized product states. Condition can then be 
written: 







(10) 



for some soluton x. The value of the distance function at the extremum is D (x). If one moves away 
from the minimum by a small displacement Sx"", the value of the distance function at the minimum 
changes by: 



dx"-dx^ 



6x^6x^ + 0{6x) 



where we have used the extremal condition (fTO!) to eliminate the term linear in the variation. The 
Hessian: 



H, 



ab 



dx"-dx'' 



(12) 



provides information about how the distance functions changes as one moves away from the extremum 
in parameter space. Specifically, the eigenvectors of the Hessian indicate the directions in parameter 
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space in the neighbourhood of the extremum in which the distance function is increasing (positive 
eigenvalue), decreasing (negative eigenvalue) or unchanging (zero eigenvalue). For the extremum to 
be a local minimum, all eigenvalues must be positive, apart from the zero eigenvalues associated with 
symmetries of the system. We now evaluate the general Hessian for the system under consideration. 
Diagonal terms of the Hessian have the form (no summation on i): 



da: 



The following off diagonal terms vanish identically: 



daodai daida^ 
All remaining terms have the form: 



dbodb 



0. 



dttidbj 



D"^ = Attibjcldf 2 {xijk...Ckdi ■ ■ ■) 



(13) 



(14) 



(15) 



Using this notation the Hessian can be written: 



H 








laobo 


7aofei 


laoco 


laoci 







laibo 


7aifei 


7a ICO 


laici 


Iboao 




no 





76oco 


76oCi 


Ibiao 


Ibiai 





n^ 


Ibico 


76ici 


Icoao 




7co6o 


Icobi 


Tco 





Iciao 


Iciai 




Ic^b^ 





n^ 



\ 



(16) 



We wish to study the unnormalized geometrical entanglement of target states that are permutation 
invariant. We make the ansatz that the closest product state has the same symmetry and restrict 
consideration to evaluating the Hessian for product states that , like the target state, are maximally 
permutational invariant: 



ai 



bo 
bi 



which gives: 



Na=N^ 



B 



«0 , 

ai , 



AT := On + ai 



(17) 



It is important to note that in general the closest product state may not necessarily have the same 
number of symmetries as the target state. It was proven in [24] that the closest normalized product 
state to a permutationally invariant target state is also permutationally invariant [31]. This has not 
been proven for unnormalized product states, so that (fT7j) is for the moment merely an ansatz. In 
order for it to provide a correct measure of the geometrical entanglement, the ansatz must not only 
provide an extremal solution to the variational equations, but as discussed above it must correspond 
to a local minimum. In section I VI I we show below that for the permutation invariant target states 
under consideration, the eigenvalues of the Hessian evaluated at extrema of the form ( ITTl) are indeed 
positive. This will prove that maximally symmetric unormalized product states do provide a local 
minimum to the distance function. It is beyond the scope of the present paper to prove that these 
are in fact global minima. 
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Equation ([6]) and the corresponding equations obtained by differentiating respect to the variables 
bj, Cfc, ■ ■ ■ form a set of 2q equations which depend on the 2q variables {oj, bj, ■ ■ ■}• When we use 
the ansatz in Eq. (fT7|) . 2(g — 1) of these equations are automatically satisfied, and the remaining two 
equations determine the values of the two paramaters and ai which correspond to an extremal 
solution. 

Using ( ITTl) . the diagonal terms in ( fT3|) are all equal: 



r,, = 2N'i-' =: r , xe{a,b,c,-- ■} 



(19) 



Each of the off-diagonal terms in f llSp is equal to one of three terms which we denote 700, 7oi or 711 
(because of the symmetry of target state (Eq. (fT7j) ). we have 710 = 7oi)- We can write all three of 
these terms collectively as (no summation on i and j): 



2 a 



a-2-{k+l+-) k+l+ 



*0 



a 



T 

T 



r 



7oo 7oi 
7io 7ii 



M 



7oi 



We define the 2x2 matrices: 
T 

and the variables: 

701^1 := - 7ii) , 701^^2 := (7oo - ^) , 7oi 
Using this notation we can write the Hessian in Eq. (fT6|) as: 



Xijk... 



V2 
1 



-Vi 



7oi 



H 



( ^ 


r 


r 


r 


T 


r 


r 


r 


T 



V 



(M + T) 
(M + T) 

V 



(M + T) 
T 

(M + T) 



(M + T) 
(M + T) 
T 



(20) 



(21) 



(22) 



J 



(23) 



IV. EIGEN-SYSTEM OF THE HESSIAN 



Substituting (fT7|) into the equations of motion (jS]) we obtain: 









Xijk-ajak ■ ■ ■ 


= 0, 




i = 


-> 


a^N'^-^ - 


Xoofe-ttottA: • ■ ■ 
■ = —{aoN^- 


- Xoifc-«i«fc ■ ■ 


= 0, 






Xoofe-«fc ■ 


- Xoifc-«i"fc ■ 


■■), 


i = 1 


-> 


aiAr«-i - 


Xiofc- aoafc ■ ■ ■ 
■ = — (aiAr<?- 


- Xiifc- «iafc ■ ■ 


= 0, 




-)■ 




- Xiofc-ttottfc ■ 


■■)• 



Similarly, substituting fll7p into f[T51) we obtain: 
7ij = 4aiajA^^"^ - 2xijfc-afca/ ■ ■ ■ , 

700 = 4aoA^^"^ - 2xoofe-afc ■ ■ ■ , 

701 = 4aoaiiV^~^ - 2xoifc -afc ■ ■ ■ ^ Xoifc-afc ■ 
711 = AalN"-^ - 2xiik-ak ■■■ ■ 



-^(7oi-4ao«iiV^-2) 



(24) 
(25) 



(26) 
(27) 
(28) 



Substituting ([2l]), ([25]) and ([27]) into ([26]) and ([28]) we obtain: 



Too 



7ii 



«o 



ai 

T 

oo 
2 



(701 - 4aoaiiV''-')] 
(4aoaiA^'"' - 7io) 



^ 7oi 

ao 



(29) 



"0 

/ 



-7io ■ 



Using these expressions and ( 122]) we obtain: 

^^2 



«0 

f^l = — , 



"0 



-1. 



(30) 



The result in (150]) allows us to obtain analytic results for the eigenvectors and eigenvalues of the 
Hessian in Eq. ( l23l) . The derivation is given in Appendix [A] and the results are listed below. The 
eigenvectors have 2q components and we write them as lists of q 2-component vectors: 



Vi = ((t;i,l),(t;i,l),(t;i,l),---)' 



V2 



ii-vi, 



(^^1,1), (0,0), 
(0,0),(t;i,l), 



f 



V^ = ((-v,,-l),{0,0),---{vul)) 

Vl = ((-^2,-l),(t^2,l),(0,0),-. 

V2 = ((-^^2,-l),(0,0),(t;2,l),-. 



eigenvalue = qr 



eigenvalue = 



(31) 



(32) 



f 



V^2,-i = ((-^;2,-l),(0,0),---(t;2,l)) 



eigenvalue = 7oi(fi — ^2) 



(33) 



eigenvalue = qr — {q ~ l)7oi(^'i — ^'2) (34) 



^2, = ((t'2,l),(t^2,l),(t'2,l),---) 

We describe the content of this equation in words. 

1. The first eigenvector is called Vi. It has eigenvalue qr. 

2. There are (g — 1) eigenvectors labelled {V2, ■ ■ -Vq} which have (— fi, —1) in the first position, 
(f 1, 1) in any one of the remaining (g — 1) positions, and (0,0) in all remaining positions. All of 
these eigenvectors have eigenvalue 0. 

3. There are (g — 1) eigenvectors labelled {V^+i ■ ■ ■ V2g-i} with the same form as the eigenvectors 

■ ■ ■ H), except that Vi is replaced by V2. They have eigenvalues 701 ("Wi — ^2). 

4. The last eigenvector is labelled V25 and has eigenvalue qr — {q — 1)701 (fi — f2). 

Using V1V2 = — 1 (see Eq. (l30l)). it is easy to see that all pairs of eigenvectors are orthogonal, 
except for pairs which have the same eigenvalue. We can construct a completely orthogonal set of 
eigenvectors using the slater determinant, in the usual way. 



V. INTERPRETATION OF EIGENVECTORS 



Now we consider the physical interpretation of the eigenvectors of the Hessian. We expect that 
these eigenvectors correspond to some kind of propagator-like normal modes. We can represent the 
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coefficients of the original product state as a six component vector: 

V^nitiai = (flo, Ol, bo, hi, Co, Ci, " " ■) = (("o, «l), («0, «l), («0, «l) " " • (35) 

We consider translating this vector by an infinitesimal amount in the direction of each of the eigen- 
vectors of the Hessian. This produces a new vector: 

K = initial + eV^^ , ie{l,2,...q}. (36) 
First we look at the eigenvector Vi. Using f l30|) we have: 

Vi = c((ao,ai), (ao,ai), (ao,ai) ■ ■ ■) , (37) 
and from ( !35|) and (136|) we obtain: 

= ((Aao, Aai), (Aao, Aai), (Aao, Aai) ■ ■ ■ )^ , A = (1 + ec) . (38) 

We conclude that translating the product state by an infinitesimal vector in the direction of Vi is 
equivalent to scaling the components of each qubit in the product state by A. 

Next we look at the eigenvectors {V2, ■ ■ ■ Vq\. Using f l30l) the V2 eigenvector can be writtenj32| 

V2 = c((-ao, -ai), (ao, "i), (0, 0) ■ ■ ■ )^ , (39) 
and from ( I35p and (136 p we obtain: 

= ((ao(l - ec), ai(l - ec)), (ao(l + ec), ai(l + ec)), (oq, ai) ■ ■ , (40) 
~ ((-^ao, -^tti), (Aao, Aai), (ao, ai) ■ ■ 0^ ^ A = e"*"' . 

We conclude that translating the product state by an infinitesimal vector in the direction of V2 
is equivalent to scaling the components of the first qubit in the product state by 1/A, and the 
components of the second qubit by A. From Eq. ([5]), the distance measure is clearly invariant 
under this transformation. It is clear that the eigenvector V3 can be treated in the same way, and 
corresponds to a scaling of the first and third qubits, and similarly for the eigenvectors V4 to Vq. 
We conclude that the zero eigenvalue corresponds to the eigenvector that points in the direction of 
a symmetry of the distance function. 

Now we consider the eigenvectors Vg+i to V2q-i- Using (15U]) V^+i can be written: 

Vq+i = c((ai, -ao), {-ai, ao), (0, 0) ■ ■ ■) , (41) 
and from f l35|) and (136|) we obtain: 

Vq^i = {R-ce{ao, ai), Rce{ao, ai), {ao, ai) ■ ■ ■)'^ , 
where Rg is the 2x2 matrix that generates a counter-clockwise rotation through an angle 6: 

^^=f"%' • (42) 

\^ sm 6 cos J ^ ' 

We conclude that translating the product state by an infinitesimal vector in the direction of Vq+i is 
equivalent to rotating the first qubit in the product state clockwise by an infinitesimal angle 6 = ec. 
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and the second qubit in the product state counter-clockwise by the same angle. It is clear that the 
eigenvectors to V2q-i can be treated in the same way. 

Finally we consider the eigenvector V2q. Using f l30|) it can be written: 

V2q = c((-ai, ao), (-«!, ao), (-"i, "o) ■ ■ ■) , (43) 
and from ( I35p and (136 p we obtain: 

which shows that translating the product state by an infinitesimal vector in the direction of V2q is 
equivalent to rotating each qubit in the product state counter-clockwise by an infinitesimal angle 
9 = ec. 

We conclude that the translation in Eq. (1360 has a simple interpretation, for each of the eigenvec- 
tors of the Hessian. 



VI. ANALYTIC SOLUTION 



In this section we consider a general class of permutation invariant target states that are consistent 
with the symmetries of the ansatz we are using, i.e. Eq. (fT7|) . We show that all eigenvalues of the 
Hessian (except the zero eigenvalues that correspond to a trivial scaling symmetry of the distance 
function) are positive, which means that the extremal solution is a local minimum, and therefore can 
be interpreted as a measure of the entanglement of the target state. In Appendix [B] we show that if 
we choose a target state that does not respect the symmetries of the original ansatz, it is not true 
that the extremum of the distance measure corresponds to a local minimum. This result shows that 
in order for the distance function to have a physical interpretation as a measure of entanglement, 
it is necessary to choose a target state that is consistent with the original ansatz, which is what we 
expect. 

We construct a target state that is symmetric under the interchange of any 2 qubits by including 
all possible permutations of p entries of "1" and q — p entries of "0". We define the normalization 
factor: 



All non-zero values of Xijk... are equal to A. The non-zero components correspond to the indices: 

ijk . . . = P( pO . Qll . 1 ) total number of terms is ( ^ j (45) 

q-p p 

For example, if g = 3 and p = 1, A = l/y/S and the permutation invariant state is: 

IV') = ^ (|0) ® |0) ® |1) + |0) ® |1) |0) + |1) ® |0) ® |0)) (46) 
v3 

The non-zero components are xooi = Xoio = Xioo = 1/ V^- 
Using ([6]) and f lTTl) the extremal solution satisfies: 
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For the simple example in (H6] )these equations are: 



73 VSai 



These are readily solved to yield N = al + al = (4/9)^/^ and af = (3/4)A^ 
In general, rearranging the equations gives: 

/ \ Y+l / \ z / . \ w 



4 



5 



where we have defined the exponents: 



.Y = 2+^, r = f^, Z = ?q. W^-^. (50) 



Using these results it is straightforward to show: 



Oil \ p ao«i V P \ Q ^ P 

These solutions can be used to obtain an analytic result for the distance measure (25| : 

Note that the above agrees with the value obtained earlier for the simple case q = 3, p = 1. 
Substituting (jS]) and 1^ into (EH) we have: 



701 = 4aoa,N'^-' - 2 af'-'a^'A ■ ( ^ _ ^ ) , (53) 



2 o Vi(l-P)P .T,-l 



q- 1 



where we have used f l49l) and f lSTj) in the last line. 
Substituting (150]) . (IFIi) and (1551) into (PT]) -( IM|) the eigenvalues are: 



ei = gr positive definite , (54) 

62 = 0, 

63 = r I 1 ) positive definite for g > 2 , 



64 = 2r positive definite . 

The result is that all eigenvalues are positive definite for q > 2 qubits, except for a set of zero 
eigenvalues that correspond to a trivial symmetry of the distance function. We conclude that the 
extremal solution corresponds to a minimum. The fact that the eigenvalues in f l54p are independent 
of p is a consequence of restricting to a target state that does not mix values of p. A target state 
constructed as a linear combination of different values of p would still satisfy the symmetry of the 
ansatz in f lT7|) . but the analytic solutions in this case are much more complicated. 
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VII. CONNECTION WITH SCHMIDT DECOMPOSITION 



In this section we will show that the Schmidt decomposition of a general unnormalized product 
state yields a very convenient basis for the evaluation of the geometrical entanglement of an arbitrary 
target state. We start with the product state in Eq. ([3]). We construct a reduced density matrix by 
tracing over all qubits except the first: 

Afi A,... = |/Vb/Vc...](^«»/°?). (55) 

We find the eigenvectors Ua^ and eigenvalues aa^ for the reduced density matrix Ma'- 



wi') = with or(i) = NaNbNc (56) 



T 



(2)=f_^,ll withap)=0 

ao 



Clearly these eigenvectors are orthogonal. We denote the normalized eigenvectors with hats and 
construct the unitary singular matrix: Aji = {ua^)j. 

We follow the same procedure for all other reduced density matrices. For example, tracing over 
all qubits but the second gives: 



= [NaNc ...]( j , (57) 



.f)=(-^,l)'withaf 



We construct the singular value decomposition (SVD) of the matrix representation of the product 
state implied by Eq. ([3]): 

cf>,jk... = [Ai.B,yCk,...]f:,y,,„, A =-^(^2 ao) ' ^ " 7^ ( 6i bo) ■■ ^^^^ 
It is easy to show that: 



E,,,... = [A-^BT^'C,^' . . .] 0,,,... = a 6u6i,6,k a := Sm... = y^N^NnNc ■■■ . (59) 

For example, for a 2 qubit state, using ([3]) and (l58l) we have: 



aibo aibi J y/N\ \ -^i ao J y/N^ \-bi bo ' 



Substituting into (l59l) we obtain: 



S = V^(J°). (61) 
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We can write the target state in the same basis: 

Xijk... [^ixBjyCfcz ■ ■ •] ^xyz... ) (62) 
^xyz... = [■^xi^yj'^zk ■ ■ ■] Xijk... ■ 

The matrix S is messy, but we don't need to use the form of S to estabhsh the connection between 
the geometric definition of entanglement in section [T] and the Schmidt decomposition. 
From Eqs. and (l62l) we have: 



D' = {tP- 01^ - 0) = (S - S) • (S - S) , (63) 

where we have defined M ■ N := Mij^^ Nij^^^^ for arbitrary tensors M and A^. Using (159|) we have: 

= l + a^ - 2(tSiii... , (64) 

where we have used = S ■ S = 1 since the target state is assumed to be normahzed. 

We can find the value of a that minimizes this distance by solving: 



Id ^2 

2 da 



. (65) 

First we show that Em... is independent of a. It is simple to see this in polar co-ordinates: 

ao = a cos 9a , ai = a sin 9a ; bo = b cos 9b , bi = b sin 9b; cq = c cos 9c, ci = c sin 9a ■■ ■ (66) 
Using the notation: ii = a , i2 = b , = c, . . . the distance function has the form (see Eq. (^): 

q q-l 

1)2 = l-2v/^a^sin(^,J J] cos{9i^) + a\ (67) 

The middle term on the rhs of this expression is a sum of terms of the form 
(sin 9a cos 9b cos 9c cos 9d- ■ ■) + (cos 9a sin 9b cos 9c cos 9d- ■■) + ■■ ■■ In each term there is 1 sine factor 
and {q—l) cosine factors, and there are q terms which correspond to q different choices for the location 
of the lone sine factor. The point is that (1^7|) depends on the g + 1 variables 9i. {j G {1, 2, ... g}) 
and cr, in contrast to ([5]), which depends on the 2q variables {ao,ai,bo,bi . . .}. The Hessian is a 
(g + 1) X (g + 1) dimensional matrix that does not have zero eigenvalues, since the symmetry which 
produces the zero eigenvalues (see section |V]) has been removed by the switch to polar co-ordinates. 
We show this explicitly in Appendix [Cl 
In polar co-ordinates the matrices in (158|) have the form: 

cos 9a — sin 6'a \ „ / cos 9b — sin 9b 



A= ^r^^" , 5 = pr^^" (68) 

\^ sm 9a cos 9a J \ sm 9b cos 9b J 

From fl62|) and f l68|) it is clear that Sm... is independent of a. Using this result flMl) and fl65|) becomes: 

1 d 



= ^ <Tc = Siii..., (69) 
which gives: 

Dl = l- al . (70) 
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To compare with the geometric distance measure we look at the angle between (f) and ^p which is 
defined in the second line of Eq. ([8]) : 

cos " = /, n ,^ /; , I ,^ ^ /- ~ ^ ~^ '^'^^ = O^c , ('1) 



and thus we have 

Dl = l- cos^ 6^ = 1- , (72) 

which agrees with the result in the first line of Eq. ([8]). 

Recall that in order to complete the calculation of the geometrical entanglement, it is necessary 
to extremize Sm with respect to the remaining parameters describing the unnormalized product 
states. Using polar coordinates these are the angles 6a, Ob--- that appear in (l^Sj) . While this will 
yield a complicated set of non-linear equations in general, in the case of the permutation symmetric 
target and product states of the previous sections things simplify considerably. In particular, for the 
product states to be permutation invariant, all the angles must be the same: 

e:=ea = eb = -- (73) 

Given Eq. (l62!) and the form (H5|) of Xijk... in the case of a permutation invariant target state, we 
have that: 

Sin... = -^(cos(^))''-P(sin(0))^ (74) 

The extremal condition dUm ^ JdO = then yields the solution: 

tan^r^,) = (75) 
q-p 

Substituting this back into Sm gives the minimum distance of 1 — with: 

^ p / \ (q—p) 



.^ = iV. = EV=i;im (l-'-l (76) 



in agreement with Eq. fl52|) . 

For clarity we illustrate the above using the example in Eqs.( l46p of the previous section, namely 
g = 3 and p = 1. In this case we have: 



A-' = B-' = C-'=i ''''' ''''' ) . (77) 

' — sm u cos a ' ^ ' 



Thus, 



Sill - ^u^ij'^ikXijk 

= cos 9 cos 9 sin 9 xooi + cos 9 sin 9 cos 9 Xoio + sin 9 cos 9 cos 9 xioo 

= 3 X ^cos^^ sin^ (78) 
v3 

as in ( FM|) . Extremizing this expression with respect to variations of 9 yields = 4/9, as expected 
from the general solution (|76ll . 
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VIII. CONCLUSIONS 

We have studied a generalisation of the usual geometric measure of entanglement of pure states 
using the distance to the nearest unnormalised product state. When the target state has a large 
degree of symmetry, one can find general examples for which one can solve the system of non-linear 
equations analytically. For these solutions, we have proven that all eigenvalues of the Hessian are 
positive, which means that the extrema are local minima. This provides a local (as opposed to global) 
version for unnormalized product states of the proof by Hubener et al [24] for normalized product 
states that permutation symmetric product states minimize the distance to permutation symmetric 
target states. In addition, we have shown that the conditions that determine the extremal solutions 
for general target states can be obtained directly by parametrizing the product states via their 
Schmidt decomposition. 

Our results verify that the distance measure we have defined in j25| is a meaningful measure 
of entanglement. In order to show that it is potentially useful for physical systems, one needs 
to analyze it in more general settings such as multi-partite target states with less symmetry and 
mixed, as opposed to pure, states. This work is currently in progress. 

Acknowledgements This paper is dedicated to the memory of our collaborator Randy Kobes, who 
passed away in September, 2010. We gratefully acknowledge the support of Natural Sciences and 
Engineering Research Council of Canada. G. Kunstatter thanks Robert James for useful discussions. 

Appendix A: Derivation of the eigenvectors and eigenvalues of the Hessian 

It is straightforward to verify the eigenvectors and eigenvalues in Eqs. fl3T]) - flM|) . The eigenvalue 
equation can be written: 

row„ ■ Kn = E„ Kn[n] , (Al) 

where row„ is the n-th row of the Hessian, Em is the m-th eigenvalue, and is the n-th component 

of the m-th eigenvector. Because of the symmetry of the Hessian and eigenvectors, we only need to 
look at the first 2 rows, and the eigenvectors Vi, V2, Vg^i and V2q. We write out the first two rows 
of the Hessian in Eq. (125]) : 

rowi = ((r, 0), {'J01V2 + r, 701), {'J01V2 + r, 701), (701^^2 + r, 701) ■ ■ ■ ) , (A2) 
row2 = ((0, r), (701, -701W1 + r), (701, -701^1 + r), (701, -701^1 + r) ■ ■ ■ ) . 

Contracting these rows with the eigenvectors Vi, V2, Vg^i and we obtain: 

rowiVi = (t;i,0)-(r,0)^ + (g-l)[(t;i,l)-(7oit;2 + r,7oi)^] (A3) 
= qrvi + (g - 1)701(1 + V1V2) = gr(fi) = Ei Vi[l] , 

row2yi = {vi, 0) ■ (0, rf + {q - l)[{vi, 1) ■ (-7oii;i + r, 701)^] = qr = Ei Vi[2] , 

rowiV2 = (r, 0) • {-Vi, -1)^ + (701^^2 + r, 701) ■ {vi, 1)^ = 7oi(l + V1V2) = = E2 V2[l] , 

row2l^2 = (0, r) ■ {-vi, -if + {--foiVi + r, 701) ■ {v^, 1)^ = = E2 V2[2] , 

rowiVg+i = (r, 0) • {-V2, -if + (701^^2 + r, 701) ■ (^2, 1)^ = 7oi(l + vj) 

= 70l{Vl - V2){-V2) = Eq+i Vq+l[l] , 
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rowaVg+i = (0,r) ■ {-V2, -l)'^ + (-701^1 + ^,7oi) ■ (^^2, 1)^ = loi{vi - t^2)(-l) = Eg+i Vg+i[2] , 

TOWiV^g = {V2,0) ■ {T,Of + {q-l)[{v2,l) ■ (701^^2 + r,7oi)^] 

= qrv2 + (g - 1)701(1 + V2) = {qr - (g - 1)701 (^^i - t'2))(f2) = ^2g V2q[l] , 

TOW2V2q = {V2,0) ■ (0, t)"^ + " 1)[(^^2, 1) " (-701^1 + ^,70l)^] 

= {qr - (g - 1)701 (^1 - ^2)) = E2q V2g[2] , 

Appendix B: Symmetry Violating Example 

Consider the target state: 

IV') = (|1100 ...0) + |0110 ...0) + |0011 ...0) + ... + |00 ... 11) + |10 ... 01)) (Bl) 

This target state does not correspond to the state in section |VT] with p = 2, because in this case Xij... 
is non-zero only when there are two adjacent I's (using periodic 'boundary conditions'): 

Xiiooo - o = Xoiioo - o = Xooiio-'-o = Xoooii---o = ■ ■ ■ = Xooooo---ii = Xioooo---oi = —p^ ■ (B2) 

The target state is not symmetric under the interchange of any two qubits, and therefore it does not 
respect the symmetry of the ansatz in Eq. (fTTI) . Using ([6]) and ( fT71) the extremal solution satisfies: 

N'^-^ = ^^^alal\ m-^ = —al-\ (B3) 

yg 

Rearranging gives: 

2 N{q - 2) 2N 

= , a^ = . (B4) 

q q 



Using f l20|) we have: 



701 = Aa^a^m-^ - ^a^ , (B5) 



"0 V g - 2 

where we used ( 1B2I) in the first line, and (1B3I) and (1B4I) in the second line. Substituting (l30l) . (1B4I) 
and ( 1B5I) into (l3T])-f lMl) the eigenvalues are: 

ei = qr positive definite , (B6) 

62 = 0, 

63 = r ( 2 — - ) positive definite for g > 3 , 

V 2(g-2)y 

64 = — r ( - — - — ^— T— I negative definite for g > 6 . 

V 2(g-2) J 

The existence of a negative eigenvalue means that the extremal solution is not a local minimum. 
This result is not unexpected, since the target state is not consistent with the symmetry required 
by the ansatz in Eq. ([1] 
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Appendix C: The Hessian in Polar co-ordinates 



Using the solution ( !76|) . defining N = A^^ and setting p = 1 for simplicity we obtain: 

q-l 



N 



9 = cos 



-1 



1 



which reproduces 



1 — N, in agreement with fl59l) and (|70 



The Hessian has the form: 



H 



B Z Z Z Z . 

Z M X X X . 

Z X M X X . 

Z X X M X . 



(CI) 



(C2) 



^ = ^^' ^ ^ 2N (C3) 

Z^^D^ - --^cosn.) + (l-.)sm^(.)cos-(.) 

M = ^Z.^ ^ 2y^^/Nsin(^)cos-^(^) ^ ^ 

X = ^ " " ''""'"'^^^ " ^ ''""^^^ cos^-'(^^)) 

^ - (^2- (g-2)(g- j N. 

In each line of ( ]C3p . the first arrow indicates the results obtained using the ansatz (173|) . and the 
second arrow indicates that we have used the equation of motion ( ICip . It is straightforward to 
calculate the eigenvalues: 

ei=B, 62 = 63 = . . . = = M - X , e^+i = M + 4X . (C4) 

As expected, the zero eigenvalues have disappeared, and we have 3 distinct eigenvalues, as before. 

It is straightforward to see how the eigenvalues in (IC4p are related to those in f l5T]) - (lM|) . We look 
at one example. Using the chain rule we have: 

D = TT^g - 1)700 + - 1)701 + ^^TijT" (C5) 



dai ( dao dai dai 

The left-hand side of this equation is the eigenvalue Ci in (1C4I) . The right-handside is a function of 
the eigenvalues in (13T]) - (1M|) . Using f l29|) and f pOj) we obtain: 
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The derivatives can be calculated directly: 



dN 2g3/2 ' 2g3/2 

It is straightforward to show that these results satisfy: 



Substituting we reproduce the first equation in flC4p . 



^11^' 1^ = 2q vi{vi - V2)tN . (C7) 
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